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DISPOSITIF DE SIMULATION DU MONDE REEL PAR TRAITEMENT 

ASYNCHRONE ET CHAOTIQUE 

5 

L'invention conceme la simulation informatique du monde reel, et notamment de F evolution 
temporelle de milieux objets de phenomenes physiques et/ou chimiques et/ou biologiques. 

10 Coimne le salt Thomme de VeaX de tres nombreux objets (animes ou inanimes) ou 
substances interagissent avec les milieux (ou media) dans lesquels ils sont places. Or, 
lorsque Thomme con9oit un nouvel objet, comme par exemple un bateau, ou bien lorsqu*il 
tente de comprendre un mecanisme reactionnel, comme par exemple le developpement d'lm 
cancer, il lui est g6n6ralement impossible de prendre en compte toutes les interactions, ou 

15 du moins les principales, notamment lorsqu'il utilise des maquettes. Par consequent, il arrive 
frequemment que Tobjet final ne procure pas les resultats escomptes, ou que Texplication 
avancee ne corresponde pas a la realite, en raison de la complexite et/ou de la combinaison 
des interactions presentes dans le monde reeL Ainsi, un bateau pent s'averer incapable de 
supporter certains etats de mer^ alors meme qu'il avait fait Tobjet d'etudes prealables 

2 0 poussees, y compris dans un bassin de carene. 

Afin de mieux anticiper le comportement final d*un objet ou produit, ou de mieux 
comprendre im mecanisme reactionnel, de nombreux industriels et de nombreux chercheurs 
utilisent des dispositifs de simulation informatique du monde reel. 

25 

Ces dispositifs de simulation reposent generalement sur une modelisation des phenomenes 
physiques et/ou chimiques et/ou biologiques sous la forme d'equatLons difFerentielles. En 

Fabsence de solution anaiytique, la simulatiorLinformatique doit s'appuyer sur des methodes 
iteratiYes. qui passmt d'un etatcourant du s3/steme- c-onsidere a un etat suivanlde ce meme 

3 0 syBttnm. CMque s\^tcni£ ai^i gCii^§rai2mmi decoBipose en sous-systknes^u mailleadoni ies 

1 ^ol jt: ?Tiz.rrcr rctr* j:: sobi cdculc-zL en pro^ik. j= :> neon r-mhrnm, Fazctrnz'^qi^zm^V ii2a 
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tions des parametres de simulation en cours de travail, ni la combinaison de modeles de 
natures difKrentes, ce qui nuit a sa generalite d' application. 

L'invention a done pow but d'ameliorer la situation. 

EUe propose a cet effet un dispositif de simulation du monde reel, propre k etre implante 

dansunordinateurcapabledesupporterenmodemulti-tachesuneprogrammationparobjet^ 
active, repr&entatifs de systemes h simuler. 

Ce dispositif se caracterise par le fait qu'il comprend un logiciel de simulation par objets de 
revolution conjoiute de certains au moins des objets actives, comportant : 

- des objets d'etat contenant chacun au moins une domiee d'espace et/ou de temps et/ou au 
moms une donn^e de propriete, definissant un 6tat courant, 

- des objets d'interaction contenant chacun la designation d'au moins I' un des objets d'^ 
et d'au moins une fonction applicable ^ I'un au moins de ces objets d'dtat, et definissal a 
chaque instant la topologie du systdme simuie, f 

-ungestionnairedesimulationcapabledetravaillerparsequencessurunese^^^^ 
d'interaction, et d'activer chaque objet d'intet^ction une unique fbis lors de cha^^ 
sequence, selon un ordre variant de fa^on au moins partielkment aieatoire d'une sequence 
h I'autre, de maniere a appliquer chacune de ses fonctions k I'etat cour^t de chaque objet 
d'etat qu'il designe pour faire evoluer son etat vers un nouvel etat courant 

En d'autres tennes, le dispositif selon l'invention fonctiom>e selonunmode asynchrone du 
fait que les etats respectifs des objets actives varient les uns apr^s les autres au sein de 
chaque sequence, compte tenu des etats respectifs des autres objets actives, et chaotique du 
fait que I'ordre de traitement de chaque objet active varie de fa9on aieatoire d'une sequence 
i I'autre. 11 est possible d'activer ou supprimer a tout moment (c'est-a-dire en temps reel) 
un ouplusieurs objets, afin de modifier les conditions de travail et/ou le systeme simule sans 
qu'al faxlle recommencer integralement la simulation, ce qui conTere au dispositif un 
veritable caractere interactif. 

Lelogicieldesimulationpeutenoutrecomprendredesobjetsd'interactionintemecontenant 
la designation d'un seul objet d'etat et d'au moins une fonction applicable a cet objet d'etat. 
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et des objets d'interaction mutuelle contenant la designation d'au moins deux objets d'etat 
et d'au moins une fonction applicable a des domiees de propriete de ces objets d'etat 
designes, 

5 Par ailleurSs certains au moins des obj ets d' etat peuvent comprendre une donnee de propriete 
qui represente une grandeur intensive, et^ou certains au moins des objets d'interaction 
peuvent comporter une fonction faisant intervenir une grandeur extensive ou intensive. 

On pent egalement utiliser des objets d'etat qui pour certains au moins d'entre eux, 
1 0 comprennent des sous-objets d'etat, ainsi qu'eventuellement des sous-objets dMnteraction 
operant sxrr ces sous-objets d^'etat 

En outre, le logiciel de simulation pent comporter des classes d'objets qui definissent des 
structures d'objets d'etat et d'objets d'interaction* Les objets d'etats et d' interaction sont 
15 alors derives de ces classes par instanciation, 

De plus, le logiciel de simulation pent comporter un ordonnanceur capable de fonctionner 
soit selon rm mode en temps reel, dans lequel il fonctionne selon une frequence choisie, soit 
selon mx mode en temps virtuel, dans lequel il fonctionne de fa^on periodique mais pendant 
2 0 des durees variables d'une periode a Pautre. 

D'autres caracteristiques et avantages de Tinvention apparaitront arexamen de la description 
detaillee ci-apres, et des dessins annexes, sur lesquels : 

- la figure 1 illustre de fafon tres schematique, sous la fomie de blocs fonctiomiels, un 
25 ordinateur equipe d'un exemple de realisation d'un dispositif de simulation selon 
rinvention, et 

■* la figure 2 est un organi gramme detaillant un exemple de fonctionnement d'un dispositif 
selon rinvendon; 
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Comme cela est illustre schematiquement et de fa9on fonctionnelle sur la figiire 1, mi 
dispositif de simulation selon rinvention D peut-8tre installe dans xin ordinateur C 
comportant un systeme d' exploitation OS et des moyens de traitement et de calcul CPU 
adaptes a un fonctionnement dans un mode multi-tdcheSj tel que celui offert par I'environne- 
5 ment oRis decrit notanunent dans le document "Systemes multi-agents'', pages 499 a 524, 
RSTI - TSI, 21/2002* Un tel environnement multi-tSches est particulierement bien adapte 
kla. programmation par objets actives, par exemple en langage C++ ou Java. L'environne- 
meht multi-tdches oRis est couple, comme le dispositif D qui est illustre, k un compilateur 
(ici nomme "compilateur de programmation objef ). Par ailleurs, T environnement oRis peut- 

10 etre couple, coixmie le dispositif D qui est illustre, a un traducteur en langage C++ (ici 
nomme "interpreteur de programmation objet") afin. d'ameliorer son efficacite par 
compilation. Get interpreteur pent meme etre adapte de maniere a constituer un compilateur 
en ligne dans lequel le code execute est un code compile en ligne et modifiable de fagon 
dynamique. Un tel environnement multi-taches, constituant une evolution d*oRis, est conmi 

15 sous le nom AR^Vi. ^ I 

Le dispositif D comprend un logiciel de simulation de revolution conjointe d'objets activ4's 
(ici nomme "simulateur general"). Plus precisement, ce logiciel de simulation (ou simulateur 
general) comporte des objets qui peuvent 6tre decomposes en deux groupes. 

20 

Un premier groupe (nonrnie "groupe d' objets d'etat") comprend des objets dits "d'etat" 
contenant chacun une ou plusieurs donnees d'espace (X,Y,Z) et/ou de temps (T) et^ou une 
ou plusieurs donnees de propriete (ou parametres, ou encore attributs), qui definissent 
ensemble un etat courant pour I'objet d'etat conceme, 

25 

Dans le cas le plus simple, dans lequel il n'existe qu'un seul objet d'etat active, I'objet 
d'interaction pent ne pas designer I'objet d'etat, celui-ci etant alors designe de fa9on 
implicite. 

3 0 Par ailleurs, dans le cas le plus simple, la donnee d'espace-temps d'un objet d'etat pent se 
r^duire au seul temps. 
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Certains au moins des objets d'etat comprennent une donnee de propriete qui represente ime 
grandeur intensive. 

Dans un premier exemple relatif a un reacteur chimique, un objet d' etat peut-dtre xin volume 
5 non localise et une donnee de propriete peut-etre une concentration d'une substance 
chimique. 

Dans un second exemple relatif a une veine sanguine, un objet d'etat peut-Stre une maille 
volumique localise dans Tespace, une donnee d'espace peut-etre la localisation geome- 
10 trique dans un repere choisi, et ime donnee de propriete peut-etre une concentration 
chimique ou la temperature a Pinterieur d'une maille. 

Un second groupe (nomme "groupe d'objets d' interaction") comprend des objets dits 
"d' interaction" (ou activites) contenant chacun la designation d'au moins Fun des objets 
1 5 d'etat et d'au moins ime fonction applicable a Vim au moins de ces objets d'etat designes* 
Le mot * Yonction" doit etre compris dans sa definition informatique et non mathematique, 
selon laquelle elle prend quelque chose (nMmporte quoi de defini) pour delivreruii resultat. 
Une fonction pent ainsi designer une methode, ou un procede, ou encore un cbmportement 
a mettre en oeuvre, 

20 

11 est important de noter qu'ime fonction pent dependre des donnees de propriete d' objets 
d'etat. 

Certains au moins des objets d" interaction comprennent une donnee de propriete qui 
2 5 represente une grandeur intensive ou extensive. 

Dans le premier exemple precite, relatif a un reacieur chimique, un objet d^ interaction peut- 
eire une reaction clumique qui manipule les donnees de propriete ^^concentrationCs)*"" de 
Tobiet d'etet '"'7oiumcl\ la reaction chimique peut-ete caracteris'ie p2r mi param^tre 
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mailles adjacentes, et le diffiiseur thermique peut-8tre caracterise par le parametre 
*'cohdtictivite thermique'*. 

Chaque objet d'interaction (ou actxvite) est done associe a au moins un objet d'etat par tme 
5 fonctioix» et un objet d'etat peut Stre associe a plusieurs objets d'interaction (ou activites). 
Une fonction ne modifie que les paranaetres du ou des objets d'etat auxquels elle est 
associee, sans modifier les parametres des autres objets d'etat et objets d'interaction. Par 

ailleurSj les caracteristiques des objets d'etat peuvent permettre de choisir la fonction 
associee aux objets d'interaction ou d'en modifier les caracteristiques. En outre, une 
1 0 fonction peut-etre modifiee de fa9on temporaire ou permanente en fonction des parametres 
d'un objet d'etat associe. 

Comine on le verra plus loin, un objet d'interaction peut constituer une interface definie,^ 
gen^ralement avec une orientation, entre deux milieirx dits "source" et **cible" (eventuelle- 1 

15 naent confbndus) constituant des objets d'etat. L'interface fiait alors le lien entre les deux ;v 

' * • ■ 
milieux, qui sont en liaison h un instant donne (notion de synchronicite forte), et modifie r 

leurs parametres respectifs (par exemple, l'interface peut etre une difference de densite entre ^ 

deux milieux qui provoque ime diffusion moleculaire a une certaine vitesse). 

2 0 Tout type d'interface peut Stre envisag^^ et notamment des interfaces de transport (par 
exemple de molecules, de chaleur, de courant, et analogues), et des interfaces de relaxation 
(par exemple de reactions chimiques ou nucleaires). Pai' consequent, xm milieu peut ne pas 
avoir d' extension spatiale, celle-ci etant induite, tout comme les eventuelles vitesses, par une 
ou plusieurs interfaces. Par ailleurs, une activite appliquee a un milieu peut modifier le 

2 5 milieu, comme par exemple des phenomenes de relaxation, ce qui correspond a une auto- 

activation interne au milieu. 

Le logiciel de simulation (ou simulateur general) comprend egalement un gestionnaire de 
simulation couple aux groupes d'objets d'etat et d'objets d'interaction et agence de mani^re 

3 0 a creer son propre sequenceur ou ordonnancetir (ou en anglais ''scheduler^') afin de travailler 

de fa9on sequentielle sur ime selection d'objets d'interaction dudit groupe d'objets 
d'interaction. Plus precisement, le gestionnaire de simulation est charge d'activer une unique 
fois lors de chaque sequence, sous le controle du sequenceur (ou scheduler) qu'il cree pour 
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roccasion, chaque objet d'interaction selectionne, selon un ordre qui varie de fa^on au 
moins parti ell ement aleatoire d'lme sequence a F autre, afin d'appliquer chacune de ses 
fonctions a Tetat courant de chaque objet d'etat qu'il designe de matiiere a faire evoluer son 
etat vers xin nouvel etat courant 

5 

En d*autres termes, conrnie cela est illustre a titre d'exemple sur rorganigramme de la figure 
2, Tutilisateur choisittout d'abord, dans une etape 10, parmi las premier et second groupes 
d'objetSj, un ou plusieurs objets d'etat et un ou plusieurs objets d'interaction qui se 
rapportent a ce ou ces objets d'6tat choisis, afin que le dispositif D simule revolution spatio- 
1 0 temporelle du systeme repr^sent^ par lesdits objets choisis, Cela "pre-active" chaque objet 
d^interaction choisi au sein du logiciel de simulation. 

Puis, dans une etape 20^ le gestionnaire de simulation initialise un compteur de sequences 
en pla9ant la valeur n du compteur a 1, et cree une liste d'objets d'etat et d'objets d'interac- 

15 tion. La premiere sequence (n=l) commence par une etape 30 dans laquelle le gestionnaire 
de simulation selectionne de faf on aleatoire, au sein. de la liste d'objets dMnteraction, Fun 
des objets dMnteraction choisis, ce qui I'active momentan^ment Chaque fonction de Tobjet 
d' interaction selectionne et active est alors appliquee k Tetat courant de chaque objet d'etat 
de la liste d' objets d'6tat avec lequel il est en relation, ce qui modifie 6ventuellement son 

2 0 etat en cours. L'activation et Tapplication (ou execution) de chaque fonction a chaque objet 
d'etat choisi constituent Tetape 40. 

Puis, dans une etape 50, le gestionnaire de simulation supprime de la liste d'objets d'interac- 
tion de la sequence en cours Fobjet d'interaction qu'il vient d'appliquer, 

25 

Enauite, dans une etape 60, le gestionnaire de simulation effectue xm test destine a 
determiner s'il reste d'autres objets d'interaction a appliquer dans la liste d' objets d.^'iaterac- 
tion de la sequence en cours* 
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eventuellement modifie par ractivation de I'objet d'interaction precedent (lequel ne peut 
plus etre utilise dans la sequence en cours). Le gestionnaire de simxdation reproduit ces 
operations (selection-activation, application et mise(s) a jour) autant de fois qu'il y a 
d'objets d'interaction selectionnes dans sa liste d'objets d'interaction, de sorte que chaque 
5 fonction de chaque objet d'interaction soit appliqu^e tme unique fois k chaque objet d'etat 
active auquel il est associe. Une fois cela termine, la premiere sequence (n=l), correspondant 
^ Pinstant T=0, est achevee, 

Lorsque la liste d'objets d*interaction de la sequence en cours est vide, le gestionnaire de 
1 0 simulation incremente d'une unite, dans une etape 70, la valeur en cours n du compteur de 
sequences. Bien entendu, cette etape 70 comporte xm test sur le nombre de sequences a 
effectuer. Si le nombre de sequences effectuees est egal au nombre maximal prevu, le 
gestionnaire de simulation met fin a la simulation. En revanche, s'il reste au moins une 
sequence a effectuer, le gestionnaire de simulation retoume a P etape 20 pour effectuer urie 
15 nouvelle sequence correspondant a im instant T+1, T+2, T+n. II reitere alors a chaque 
nouvelle sequence les operations pr^citees. 

La duree de la simulation, et done le nombre maximal de sequences effectuees par le 
gestionnaire de simulation, depend de Tapplication concem^e, ou du param^trage choisi par 

2 0 Tutilisateur compte tenu de Tapplication. Mais, la simulation peut gtre interrompue a chaque 

instant par Tutilisateur a Taide d'une instruction d'arrSt transmise au logiciel de simulation 
grace a une interface homme/machine de Tordinateur C. 11 est important de noter qu'une 
simulation interrompue a la demande d*un utilisatejxr peut-Stre reprise ulterieurement. 

25 En raison de ce mode de fonctionnement sequentiel et chaotique du logiciel de simulation, 
Tutilisateur peut a chaque instant intervenir dans une simulation, soit sous la forme d^un 
"avatar"' pour interagir lui-mSme avec Tobjet de la simulation, par exemple le milieu, soit 
pour aj outer a, ou supprimer de, sa selection un ou plusieurs objets d'etat et^ou un ou 
plusieiars objets d'interaction. L'utilisateur peut egalement decider de modifier au moins 

3 0 partiellement la definition (ou structure) d'un ou plusieurs objets d'etat ou d'interaction. 

Cela confere au logiciel de simulation une grande interactivity. 
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Trois exemples permettant de mieux comprendre ce que representent les objets d'etat et 
d'interactions sont donnes ci-apres. Pour eviter d'alourdir inutilement la description, les 
exemples sont simplifies et reduits a ce qui est necessaire pour permettre a rhonime du 
metier de comprendre comment on peut mettre en oeuvre Tinvention. 

5 

Un premier exemple conceme la simulation de la cinetique chimique au sein d'un r6acteur 
chimique. Le r6acteiir chimique constitue un obj et d' etat representant un milieu comportant 
des substances chimiques, et au sein duquel peuvent survenir N reactions chimiques 
constituant autant d'objets d'interaction* 

10 

L'objet d'etat reacteur chimique est par exemple associe k des donnees de propriete 
representant les concentrations (C1(T)5 C2(TX Cn(T)) a I'instant T des n substances 
chimiques presentes initialement dans le reacteur. Ces concentrations definissent ici I'etat 
(a rinstant T) de Fobjet d'etat reacteur chimique, dont on souhaite simuler revolution 
15 temporelle. 

Chaque objet d'interaction reaction chimique est par exemple caract6rise par des donnees 
de propriete representant, d*une premiere part, la vitesse de la reaction V = f(Cl, C2, Cn), 
d'une deuxieme part, les reactifs concemes par la reaction (Rl , Rn), d*une troisieme 

20 part, les produits resultants de la reaction (PI, P2, Pm), et d'une quatrieme part les 
coejffiicients stoechiometriques definissant les proportions initiales des reactifs a T=0. En 
outre, chaque objet d'interaction reaction chimique est par exemple associe a une fonction 
representant le comportement de la reaction : dans un premier temps on lit les concentrations 
(C1(T), C2(T% ,..5 Cn(T)) en cours des reactifs presents a Finstant T dans le reacteur 

25 chimique. dans un deuxieme temps on calcule la vitesse de reaction V qui depend des 
concentrations des reactifs presents dans le reacteur chimique, et dans un troisieme temps 
on modifie les concentrations des reactifs et des produits apres un temps dT a la vitesse de 
reaction^/. 
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Dans ce premier exemple, I'utilisateiir va done choisir son objet d'etat (c'est-a-dire la 
composition initiale de son reacteur chimique) at ses N objets d'interaction (c'est-a-dire les 
N reactions chimiques qui sont impliquees dn fait de la composition choisie). Le 
gestionnaire de simulation pent alors debater son traitement sequentiel. 

5 • ' 

Comme indiqne precedeimnent, il debute une premiere sequence en selectionnant de fa9on 
aleatoire Tun des N objets d'interaction choisis. Get objet d'interaction est active afin qu'il 
applique sa fonction (comportement de la reaction) a I'^tat courant de Tunique objet d'etat, 
defini ici par les concentrations (C1(T), C2(T), Cn(T)) en cours des reactifs (lors de la 

10 premiere sequence T=0). Le gestionnaire met a jour les concentrations (C1(T), C2(T), 

Cn(T)) des reactifs. Puis, il selectionne de fafon aleatoire Tun des N-1 objets d'interaction 
restants, afin de Factiver et d'appliquer sa fonction au nouvel etat courant de Tunique objet 
d'etat, II reproduit ces operations N fois de sorte que chaque fonction de chaque objet 
d' interaction soit appliquee une unique fois 4 Vobjet d^^tat La premiere sequence est alors 

15 terminee. 

Le gestionnaire debute alors une nouvelle sequence, correspondant a Tinstant T+l, si celk 
s'avdre n^cessaire, en reiterant les operations pr^citees. La simulation prend fin lorsqu*il ne 
reste plus de reactifs pour produire des produits, ou bien lorsque Tutilisateur adresse une 
2 0 instruction d'arret au logiciel de simulation, 

Un deuxieme exemple conceme la simulation de la diffusion moleculaire de M diflftiseurs 
au sein d'un espace decoupe en K mailles, 

25 Chacune des K mailles constitue un objet d^etat caracterise^ par exemple, par une donnee 
deposition (ou localisation topographique ou geographique dans un repere 3D) et par des 
donnees de propriete representant les concentrations (CI (T), C2(T), Cn(T)) a Pinstant 
T de n substances chimiques (SI, S2, Sn) presentes dans lamaille. Ces concentrations 
definissent ici Petat (k Finstant T) de la maille. 

30 

Chacun des M diffuseurs constitue un objet d'interaction caracterise, par exemple, par des 
donnees de propriete (ouparametres)repr&entant, d 'une premiere part, la substance difius^e 
Si, d'une deuxidme part, la vitesse de difiusion tridimensionnelle V = (Vx, Vy, Vz) de la 
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substance Si, et d'une troisieme part, les deux mailles entre lesquelies s'effectue la diffusion. 
En outre, chaque diffiiseur (objet d' interaction) est par example associ^ a une fonction 
representant son comportement : dans un premier temps on lit la concentration Ci(T) de la 
substance Si k diffuser dans chacinie des deux mailles^ dans im deuxieme temps on calcule 
5 les quantity a difiEuser (par exemple en utilisant la loi de Fick generalisee), et dans un 
troisieme temps on modifie la concentration du diffuseur Si dans chacxme des deux mailles 
apres un temps dT k la vitesse de diffusion V- Les donn^es de propriete representent ici les 
attributs (ou parametres) de la fonction (comportement du diffuseur). 

10 Ici chaque objet d'interaction agit sur deux objets d*etat. de sorte qu'il constitue rni objet 
d'interaction mutuelle. Par ailleurs^ chaque objet d'interactionpresente ici une structure (ou 
definition) classique reposant sur au moins une fonction designant au moins un objet d'etat, 
ainsi qu'eventuellement des parametres ou regle(sX ou encore loi(s), lies a la fonction et 
constituant des donnees de propriete* Mais, leur structure pent Stre plus complexe, par 

15 exemple lorsqu'elle comporte un ou plusieurs sous-objets d'interaction. En outre, les mailles 
sont ici de type tridimensionnel (3D), mais dans certaines applications elles peuvent 8tre de 
type bidimensionnel (2D), voire mdme xmidimensionnel (ID). 

Dans ce deuxieme exemple, Tutilisateur va done choisir ses K objets d'6tat (c'est-i-dire les 
20 K mailles de diffusion) et ses M objets d'interaction (c'est-a-dire les M diffiiseurs). Le 
gestionnaire de simulation pent alors debuter son traitement sequentiel. 

11 debute une premiere sequence en selectionnant de fa9on aleatoire Tun des M objets 
d'interaction choisis. Cat objet d'interaction est active afin qu'il applique sa fonction 
25 (comportement du diffuseur), compte tenu de ses propres attributs, a Petat courant des K 
objets d'etat (mailles), deJSni ici par les concentrations (C1(T). C2(T), Cn(T)) en cours 
des substances chimiques (SL S2. Sn) dans les differentes mailles (lors de la premiere 
sequence T=0). Le gestiormaiis met a jour les concentrations. (C I (T), C2(T), Cn(T)) des 
sutetanaes {3L 22, .... Sn) dmn cliacuns des il in2iil2£^* Puis,, il i^^leGtioim^ diLiagon 
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Le gestionnaire debute alors une nouvelle sequence, correspondant a rinstant T+l, si cela 
s'avere necessaire, en r^iterant les operations precitees. La simulation prend fin lorsque les 
concentrations (C1(T), C2(T), Cn(T)) des substances chimiques (SI, S2, .... Sn) sont 
respectivement identiques dans chacune des K mailles, ou bien lorsque rutilisateur adresse 
5 une instruction d'arrSt au logiciel de simulation. 

Dans un troisieme exemple, on peut combiner les deux exemples precedents de maniere a 
simuler revolution de la cinetique chimique et de la diffusion moleculaire au sein d'lme 
mSmeveine. 

10 

Par exemple, on choisit K objets d'etat constituant N reacteurs chimiques localises dans 
r espace 3D. Chaque reacteur chimique est par exemple associe a une donn^e de position (ou 
localisation geographique dans un repere 3D) et a des domiees de propriete representant les 
concentrations (C1(T), C2(T), Cn(T)) a Pinstant T de n substances chimiques (Si|S2, 
1 5 Sn) presentes dans la maille, Ces concentrations definissent ici Tetat (a Tinstant T) de la 

, 'V* 

mailte. 

Par ailleurs, on choisit par exemple N+M objets d'interaction constituant N reactions 
chimiques et M difiuseurs. 

20 . 

Chaque reaction chimique est par exemple associee a des donnees de propriete representant, 
d'lme premiere part, la vitesse de la reaction V = f(C 1 , C2, , , Cn), d' une deuxieme part, les 
reactifs concemes par la reaction (Rl, R2, Rn), a T=0, d'une troisieme part, les produits 
resuhants de la reaction (PI, P2, Pm), et d'une quatrieme part les coefficients 

2 5 stoechiometriques definissant les proportions initiales des reactifs a T=0. Chaque reaction 
chimique est egalement, par exemple, associe a une fonction representant le comportement 
de la reaction : dans un premier temps on lit les concentrations (C1(T), C2(T), Cn(T)) en 
cours des reactifs presents a Pinstant T dans le reacteur chimique, dans un deuxieme temps 
on calcule la vitesse de reaction V qui depend des concentrations des reactifs presents dans 

30 le reactexir chimique, et dans un troisieme temps on modifie les concentrations des reactifs 
. et des produits apres un temps dT a la vitesse de reaction V. 
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Chaque diffiisevir est par exemple associe k des doimees de propriae representant, d'une 
premiere part, la substance difftisee Si, d'une deiixifeme part, la vitesse de diffusion 
tridimensiomelle V = (Vx, Vy, Vz) de la substance Si, et d'une troisieme part, les deux 
mailles entrelesquelless'effectueladiffiision.Chaquedifiuseur est 6gal^ent,par exemple, 

5 associ^ a une fonction representant son comportement : dans un premier temps on lit la 
concentration Ci(T) de la substance Si a diffuser daas chacune des deux mailles, dans un 
deuxieme temps on calcule les quantites a difiuser (par exemple en utilisant la loi de Fick 
g^neralisee), et dans un troisidme temps on raodifie la concentration du dififuseur Si dans 
chacune des deux mailles apres un temps dT a la vitesse de division V. Les donnees de 
10 propriety repiesentent ici les attributs (ou parametres) de la fonction (comportement du 
diffuseur). 

Le fonctionnement du dispositif est alors identique k celui presente dans les exemples 
pr^cMents, mais cette fois les concentrations des differentes substances ^oluent non 
15 seulement en fonction des reactions chimiques. mais egalement en fonction des mailles oix 
elles s'effectuent. Par consequent, chaque sequence comporte N+M tirages aleatoires 
permettant d'appHquer successivement les fonctions des N+M objets d'interaction aux K 
r6actevirs chimiques. 

2 0 Bien entendu, invention pent s'appliquer k des applications beaucoup plus complexes que 
celles presentees ci-dessus a titre d'exemple iUustratif. EUe s'appUque egalement k des 
applications plus simples dans lesquelles le logiciel n'intervient que sur un unique objet 
d'etat active, a Taide d'un unique objet d'interaction interne. Par aiUeurs, dans certaines 
applications I'un au moins des objets d'etat pent presenter une structure (ou definition) 

25 complese reposant sur un ou piusieurs sous-objets d'etat (presentant la structure classique 
d'un objet d'etat), 6ventuellement associ6 a un ou piusieurs sous-objecs d'interaction 
(presentant la structure classique d'un objet d'interaction). 



1 est mamtmontfait reference 2U.iableau de I'aime-e XI. pour.uns description, au niveau. 
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10 



) 



Le code n'apparalt que dans la seconde colonne du tableau, dont la premiere colonne ne 
comporte que des reperes alphanimieriques ne faisant pas partie du code. 

L'annexe XI commence, en A, par des declarations "#mclude" qu'il est inutile de detaiUer 
du fait qu'elles sont bien connues de I'homme du m6tier. Ensuite, en B, est d^clar^ uhe 
classe interface, associ^e k la classe milieu trait^e en detail en C. Apres l'en-t6te de classe 
CI. la classe milieu comprend la declaration de m^tfaodes specifiques C20 k C30, puis en 
C4 une fonction (ou methode) "activate-, et eniin des variables protegees enoncee's en C5. 
Dans l'annexe XI, les declarations de m^odes (ou fonctions) sont allegees des defimtions 
de lype/parametres, puisque I'homme de metier retrouvera celles-ci dans I'^nonce detaiUe 
de chaque methode. 



Parmi les variables protegees indiqu6es en C5, on pent noter : 

- des "Jnterfaces", qui sont regroupees au sein d'un tableau d'un type particulier 4 
15 StlVector (classe standard C-H-), et .- % 

. ;- une variable qui est de type ArRef (classe propre ^ I'environnem 

Vientensuite le detail de laclasse interface (par souci d'homogeneite elle reste designeepaf 
B). Elle comprend une methode de definition et une methode de construction, suivies de 
0 declaration de m^thodes specifiques B20 a B31 et de la methode activate B4. Enfin sont 
definies des variables protegees B5 ou Ton note a nouveau en particulier : 

- une "jsource" qui est de type milieu, et 

- une "jarget" qui est egalement de type milieu. 

5 Les definitions des methodes de la classe milieu interviennent ensuite (en C 1 1 -01 4 et C20- 

C30)engranddetail.LesrubriquesCllaC13definissentleconstructeurde laclasse 
En CI 1 sont initialisees les variables protegees de la classe milieu. 



En CI2 apparaJt une d6finition d'une forme 3D dite "Vrml". Le code Vrml permet de faire 
de la visualisation 3D en ligne. Cette forme est definie par une chatne de caracteres 

confonnementa]asyntaxeVnnldecriteparexemplearadressehttp://www.vnnl.org.Cette 
forme est ensuite visualisee et mise en couleurs. 
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15 

La rubrique CI 3 associe au Tnilieu une activite destinee a declencher la methode activate de 
larubrique CI 4. 

La rubrique €14 est un destructe\ir de I'objet d'etat milieu. 

Les rubriques C20 k C30 incluse doraient le d^ail de la construction des methodes 
specifiques deja evoqu6es sous lamSme identite pour la classe milieu. 

La methode activate a appliquer au milieu est definie en C4. Elle permet de relier une 
1 0 indication de couleur a I'indication de densite, et plus pr^cisement de changer la couleur par 
tirage aleatoire. En d'autres termes, on modifie les attributs rouge, vert, bleu du milieu 
chaque fois qa^ activate est appliquee (ou activee). 

Interviennent ensuite des definitions pour la classe interface, avec en Bl 0-B 1 1 la definition 
15 du constructeur. Plus precisement, en Bl 0 on initialise les variables prot^g^es de la classe 
interface, tandis qu'en Bll on associe k I'interface une activite destinee k appHquer la 

methode activate definie en B4. 



20 



La rubrique B12 est un destructeur de I'objet d'interaction interface. 
Viennent ensuite les details des methodes specifiques B20 a B31 deja citees. 



La methode activate de V interface est definie en B4. II s'agit de calculer une diffusivite en 
fonction d'une densite de milieu source et d'une densite de milieu cible, afin de propager 
2 5 I'ecart entre les deux densites, si du mains cet ecart est significatif. L'ecart note "delta" est 
calcula^conmie rapplication d'une fonction difilisivite k la densite de mUieu source rhos, 
diminuee de la densite de ?Hz7iez/ cible rhoi. Pour ce faire, on demande d'abord au milieu 
source et au rAiliett cibh queUes sent leura densitss respectives. Puis,_on oalcuie la valeur 
:.bGoh3s des. diSsrences de couleurc.. Zi carte vdeur afcsolue est superieure a w^i: Taieur 
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Les rubriques B60-B63 pennettent de defmir le type des fenetres de visualisation 3D qui 
seront utilisees par rapplication, 

eest alors qu'intervient uu element important pour I'invention qui est I'ordonnanceur ou 
sequenceur (ou en anglais "scheduler^). Les rubriques S 1 0, S20, S21, S22, S30, S3 1 et S41 
pennettent d'initialiser la simulation. 

En I'esp^ce, la rubrique SIO initialise un scheduler en temps virtue!, c'est-a-dire que son 
fonctionnement n'est pas astreint a respecter le temps reel, mais chacune de ses iterations 
repr^sente logiquement, et non physiquement, une duree de une milliseconde (1 ms). Bien 
^^*^;i^l^^^''^^«^^'-peutegalementfonctiomierent^^^^ 
Iterations dure physiquement une periods choisie. 

La rubrique S20 permet d'initialiser les dimensions du milieu grtce aux arguments de a 

15 1^8«edecommandeXarubriqueS21pennetdecreeretd'initialiserlesmilieux:Larubriqtf^ 
S22 permet de creeret d'initialiser les interfaces. ^ 

Les rubriques S31 etS32r^alisentunaffichageauxiliaire,qui sefaiticiparlavoiespeciak 
de communication avec I'ecran dit "flot d'erreur". 



10 
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La rubrique S41 affiche la scene a traiter, dans des conditions de 



perspective choisies. 



5 



Enfm, la rubrique S50 est le programme principal main qui initialise le fonctionnement du 
syst^me. Plus pr^cisement, ce programme cree tout d'abord un syst^me 3D, puis il 1' appelle 
pour qu'il active la mdthode domiee a la rubrique SIO afin qu'elle cree le scheduler (ou 
sequenceur), les milieux et I'interface. 

On comprendra que la stmcture de classes, presentee ci-avant, pent etre resumee de la 
mamere indiquee schematiquement dans Tannexe X2. 

La transition entre la representation de I'amaexe X2 et le code ddtaille de Pannexe XI est 
considereecommeaccessibleal'hommedum6tier,d^lorsqu'unexempIeenaetedom.e. 
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Lnnventiontrouvedetresnombreuses applications dansdenombreuxdomain^ 
et notamment dans les domaines de la chimie, de la phaxmacie, de la physique, de 
I'aeronautique, de 1' architecture, en particulier navale (par exempie pour 1' 6tude comporte- 
mental d'un bateau ou d'une plate-forme offshore, en remplacement et/ou en complement 
des basshxs de car^e), de lamedecine, notamment dans le cadre de I'etude du developpe- 
ment et du traitement de certaines maladies (par exempie les cancers) ou de m^canismes 
r6actionnels (par exempie I'activation de I'insuline), de I'ergonomie, notaimnent pour la 
realisation d'6quipements sp^cifiquement adaptes a des personnes handicap6es, et dans le 
domaine de la circulation routifere. 

Uiavention ne se limite pas aux modes de realisation de dispositif de simulation d^crits ci- 
avant, seulement a titre d'exemple, mais elle englobe toutes les variantes que pourm 
envisager I'homme de Tart dans le cadre des revendications ci-apres. 
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Annexe XI 
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B20 


setVelocityO 


B21 


getVelocityO 


B22 


setDiffusivityO 


B23 


getDiffusivityO 


B24 ' 


setSourceO 


B25 


getSourceO 


B26 


accessSourceO 


B27 


setTargetO 


B28 


getTargetO 


B29 


accessTargetO 


B30 


1 getOtherEndO 


B31 


1 accessOfherEndO 
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Cll 



^^tivaterArRef<Activity> act, double dt); 



B5 I protected: 

ArRef<Miiieu> _soxirce; 
ArRef<Milieii> jtarget; 
double jvelocity, _diflfusivity; 

}; 



AR_Cl. ASSJ)EF(Milieu,Object33>) 

Milieu::Milieu(ArCW & arCW) 
: Object3D(arCW), 
_interfacesO, 

Jt(newJVrmlMaterialInteractorO), 

_red(1.0), 

_green(1.0), 

Jblue(l.O), 

_density(0)5 
jchangeDensity(true), 
isSovtrc€(false) 
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C12 


{ 

ArRef<VrmlShape3D> sh=new_VrmlShape3D0; 
sh->parseString(" \ 
Shape { \ 

• ' appearance Appearance { \ 
material DBF MAT Material { \ 
difftiseColorl.OLOl.O \ 
transparency 0.5 \ 

} \ 

geometry Box { \ 
size 0.9 0.9 0.9 \ 

}"); 




C13 


ArRef<Activxty> act=new_Activity(0. 1); 
act->setBehaviorCthisRefO,&Milieu::activate); 

> i 




C14 


Milie«:;-Milieu(void) 

{ 

> 




C20 


void 

Milieu::setDensity(double density) 

{ ' 
^density = density; 
_changeDensity = true; 
} 




C21 


void 

Milieu::getDensity(double & density) 
{ 

density = ^density; 

} 
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C22 


void 

Milieu::setColor(double red, 
double green, 
double blue) 

{ 

jred=Ted; 

green^green; 
_blue=blue; 

Jt->vvriteDifEuseColorField(j'ed,_green,_blue); 

} . • 




C23 


void 

Milieu: :getCoIor(double & redOut, 
double & greenOut, 
double Sc blueOut) const 

{ 

redOut= red; 

greenOut=_green; 

blueOut=Jblue; 

} 




C24 


void 

Milieu::addSource(void) 
{ 

_isSource = true; 
} 




C25 


void 

Milieu::removeSource(void) 

/ 
I 

JsSource = false; 
} 




C26 


void 

Milieu::addlnterfaG.e(ArReMnterfaGe> interface) 
\ 
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C27 


void 

Milieu: :removeInterface(ArRef<Interface> interface) 
{ 

for(unsigned int i=_interfaces.si2eO;i— ;) 

{ 

if(_interfaces[i]==mterface) 
{ 

StlVectorFastErase(_interfeces,i); 

break; 

} 

} 

} 












C28 


unsigned int 

Milieu::getNfbInterfaces(void) const 

{ 

retum(_interfaces.si2:e0); 
} 




C29 


ArConstRef<Interface> 

Milieu: :getlnterface(unsigned int index) const 

{ 

retuniC_interfaces[index]); 

) 








C30 


ArRef<Jnterface> 

Milieu::accesslnterface(unsigned int index) 

{ 

retum(_interfeces[index]); 
} 
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C4 I bool 

Milieu::activate(ArRef<Activity> act, 

double /*dt*/) 

{ 

ifCisSource) { _density = 1,0; } 
if(_changeDensity) 

{ 

setCoior(1.0,1.0-_density,1.0-_density); 

if(_density < 0.1) { setColor(l,l,l); } 
else {setColorO, 0,0);} 

*/ 

_changeDensity = false; 
} 

retum(true); 

> 



BIO AR_CLASS__DEF(Interface,ArObject) 



lnterface::Interface(ATCW & arCW) 

; ATObject(arCW), 

jsoiirceO, 

JtargetO, 
__yelocity(0), 

_diffusivity(0) 



Bll 



ArRef<Activity> actr=new_Activity(0.1); 
act->setBehavior(thisRefOAInterface::activate); 

} 



B12 



lnterface::-Iiiterface(vQid) 
{ 



B20 



void 

Int Qmceiise-f Vel ociic^<do?jhleL_Yelocm' ) 
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B21 


void 

Interface::getVelocity(doxible & vetocity) 
{ 

velocity =_velocity; 

} . 




B22 


void 

Interface: :setDififiisivity(double diffusivity) 
{ 

_difiEusivity = difiusivity; 
} 




B23 


void 

Ihterface:;ge1pijSusivity(doubIe & difiusivily) . 
{ 

difftjsivity = _difiusivity; - 
} 




B24 


void 

Interface::setSource(ArRef<Milieu> source) 
{ 

if(_soiirce) 

{ 

_source->rembveInterfece(thisRefl;)); 
} 

_source=source; 
ifl[_source) 
{ 

_source->addInterface(thisRefl3); 
} 

} 




B25 


ArConstRef<Milieu> 
Interface: :getSource(void) const 
{ 

retuTn(_soiirce); 

} 




B26 


ArRef<Milieu> 

Interface: :accessSource(void) 

{ 

retum(_source); 

} 
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B4 


bool 

Interface:;activate(ArRef<Activity> /* act */, 
double /*dt*/) 

{ 

double rhos,rhot; 
_source->getDensity(rhos); 
_target->getDensity(rhot); 
bool change^^false; 
if(fabs(rhot-.rhos)>0.01) 
{ 

double delta~_diffusivity*(rhos-rhot); 

rhos -= delta; 

rhot delta; 

change=tme; 

} 

if(change) 

_source->setDenslty(rhos); 

__target->selDensity(rhot); 

} 

retum(trae); 

} 


1 


B60 


class MyViewer : public BxaminerSD 
{ 

AR_CLASS(MyViewer) 
AR__CONSTRUCTOR(MyViewer) 
protected: 

virtual void _onKeyPress(const Viewer3D::KeyPressEvent & evt); 

}; 
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AR_CLASS_DEF(MyViewer,Examiner3D) 

MyViewer::MyViewer(ArCW & arCW) 

: Examiner3D(arCW) 

{ 

} 1 
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B62 


MyViewer: :~My Viewer(void) 
{ 

} 
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void 

MyViewer: :_onKeyPress(const Viewer3D::KeyPressEvent & evt) 
{ 

Examiiier3D::__ojQKeyPress(evt); 
if((evt,key="+")li(evt.key="-")) 
{ 

double dist^yawjpitch; 

computeCiu:sorDirection(evt.xMoiisejevt.yMouse,yaw,pitch); 
ArRef<Base3D> ray=iiewJ3ase3D0; 
ray->setLocation(thisRefO); 

ray->yaw(yaw); 
ray->pitch(pitch); 
ArRef<Object3D> found; 

getSceneO->iirstRayIntersectiori(ray,3,tme,Miiieu: :CLASS09 
foundjdistjthisRefO); 

ifl[fotind) 
{ 

ArRef<Milieu>m(static_jcast<Milieu*>(foundx_ptrO)); 

cerr « "Found " « m->getAliasO « " " « evt.key « endl; 

if(evt.key=="+'0 { m->addSourceO; } 

if(evt.key=="-") { m->removeSource(); } 

} 

} 

} 




SIO 


ArRef<Scheduler> 

simulationlDit(void) 

{ 

ArRef<Scheduler> sched==new_VirtuaiIiineScheduIer(l e-3); 
unsigned m nbX=40,nb3£=20,iibZr^2: 
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S20 I if(ArSystem::getConunandLme0.size()>l) 
{ 



strToUint(nbX,ArSystem::getConimandLineOri])- 
} 

ifi:ArSysteni::getCommandLine0.sizeO>2) 
{ 

strToUint(iibY,ArSystem::getConiinandLineO[2]); 

if(ArSystem::getCommandLine0.sizeO>3) 
{ 

strToUint(nb2,ArSystem::getCoimnandLmeO[3]); 



S21 j unsigned int x,y,z; 

StlVector<StIVector<StlVector<ArRef<Milieu> > > : 
for(x=0;x<nbX;x++) 

{ 

StlVector<StIVector<ArRef<MiIieu> > > ex; 
t.push_back(ex); 
for(y=0;y<iibY;y-H-) 

{ . 

StlVector<ArRefcMilieu> > ey; 
t.back().push_back(ey); 
for(z=0;2<nb2;z++) 
{ 

ArRef<Milieu> m=new_MilieuO; 
tback0.back0.push_back(m); 
m->setPosition(x,y,z); 
} 

} 

} 
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for(x-0;x<nbX;x++) 
{ 

for(r=0;y<nbY;y-H-) 
{ 

for(z=0;z<nbZ;z-H-) 
{ 

ArRef<Interface> i; 
double diflusivity = 0.5; 
//double velocity =1.0; 
if(x>0) 
{ 

i=new_InterfaceO; 
i-.>setSource(t[x][y][zl); 
i->setTarget(t [X- 1 ] [y] [z]); 

i-.>setDiffiisivity(diffiisivity); 

i->setVelocity(0.0); 

} 

if(y>0) 
{ 

i=new_InterfaceO; 
i->setSource(t[x] [y][z]); 
i->setTarget(t[x]Iy-l][z]); 

i->selDiffusivity(diffusivity); 
i->setVelocity(0.0); 

} 

if(z>0) 

{ 

i=new_Interface05 
i->setSource(tM[y]lz]); 
i->setTarget(t[x]Iy][2-l]); 

i->setDiffasivity(difiusivily/l 000); 
i->s€tVelocity(0,0); 

} 
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S32 


cerr << ArClass::findC'Miyeu")'">getNbInstances(tme) 
<< " milieux" « endi; 




S41 


ArRef<Viewer3D> v=NEW_MyViewerO; 

ArRef<BoundingBox3D> bbox=new_BoimdingBox3D0; 

v->getSceneO->computeBoiindingBox(bbox); 

v->setFarDista3Qce(5.0*bbox->getMaxSize()); 

v->setNearDistance(0.00 1 *v->getFarDistanceQ); 

irasigned int w,h; 

Viewer3D::getScreenSi2e(w,h); 

v->setWindowGeometry(w-400,0,400,400); 

v->setMapped(trae); 

retum(sched),' 
> 
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int 

main(int argc, 
char ** argv) 

{ 

Ar$ystem3D init(argc,argv); 
retiuii(init.sim\ilate(&siniulatioii[mt)); 

} 
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Annexe X2 
Representation simplifiee 



Classe Interface 


Constructeiir d'objets 


M6thodes specifiques 


setVeiocityO 




getVelocityO 




setDiffusivityO 




getDifftisivityO 




setSourceO 




getSourceO 




accessSourceO 




setTargetO 




getTargetO 




accessTargetO 




getOtherEndO 




accessOtherEndO 


Methode activate 


activateO 



Classe Milieu 


Constructeur d'objets 


M^thodes specifiques 


setDensityO 




getDensityO 




setColorO 




getColorO 




addSourceO 




removeSourceO 




addlnterfaceO 




removelnterfaceO 




getlSfblnterfaaeO 




getlnieiiaceO 




accessintciiaceO 
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REVENDICATIONS 

Dispositif (D) de simulation du monde reeUpropre a etre implante dans un ordinateur 
(C) agence pour supporter en mode multi-taches une programmation par objets 
actives, representatifs de systemes h simulerj, caracterise en ce qu'il comprend un 
logiciel de simulation par objets de revolution conjointe de certains aumoins desdits 
objets actives, comportant : 

- des objets d^'etat contenant chacun au nioins une donnee d'espace et/ou de temps 
et/ou au moins une doimee de propriete, definissant im etat courant, 

- des objets d'interaction contenant chacun la designation d'au moins Tun desdits 
objets d'etat et d'au moins ime fonctioii applicable a Pun au moins de ces objets d'etat, 
et definissant 4 chaque instant la topologie du systfeme simule, 

- un gestionnaire de simulation capable de travailler par sequences sur une selection 
d'objets d' interaction, et a activer chaque objet dMnteraction une imique fois lors de 
chaque sequence, selon un ordre variant de fa$on au moins partiellement aleatoire 
d'lme sequence a I'autre, de maniere a appliquer chacune de ses fonctions a Tetat 
courant de chaque objet d*6tat qu'il dfeigne pour faire evoluer son etat vers un nouvel 
etatcourant. 

Dispositif selon la revendication 1^ caracterise en ce que ledit logiciel de simulation 
comprend des objets d'interaction interne, aptes a contenir chacun la designation 
d'un seul objet d'etat et d'au moins une fonction applicable a cet objet d^etat, et des 
objets d'interaction mutuelle, aptes a contenir chacun la designation d 'au moins deux 
objets d'etat et d'au moins ime fonction applicable a des donnees de propriete de ces 
objets d'etat designes. 

Dispositif selon Tune des revendications 1 et 2, caracterise en ce que ledit logiciel 
de simulation est agence pour modifier certaines au moins desdites fonctions en 
fonction d'au moins une donnee de propriete d'au moins un objet d'etat associe. 
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Dispositif selon Tune des revendications 1 a 3, caracterise en ce que ledit logiciel de 
simulation est agence pour selectionner certaines au moins desdites fonctions en 
fonction d^'au moins une donnee de propriete d'au moins un objet d'etat associe. 

Dispositif selon Tune des revendications 1 a 4^ caracteris^ en ce que certains au 
moins des objets d'etat compreiment une donnee de propriete representant une 
grandeur intensive. 

Dispositif selon Tune des revendications 1 a 5, caracterise en ce que certains au 
moins des objets d'interaction ont une fonction faisant intervenir une grandeur 
extensive ou intensive. 

Dispositif selon Tune des revendications precedentes, caracterise en ce que certains 
au moins des objets d^'etat comprennent des sous-objets d*etat. 

Dispositif selon larevendication 7> caracterise en ce que certains au moins des objets 
d'etat comprennent des sous-objets d'interaction operant sur lesdits sous-objets 
d'etat 

Dispositif selon Tune des revendications prec^dentes, caracterise en ce que ledit 
logiciel de simulation comporte des classes d'objets definissaht des structures 
d'objets d'etat et d'objets d'interaction, lesdits objets d'etats et d'interaction etant 
derives de ces classes par instanciation. 

Dispositif selon Tune des revendications precedentes, caracterise en ce que ledit 
logiciel de simulation comporte un ordonnanceurpropre a fonctionner selon Tun de 
deuxmodes choisis parmi unmode en temps reel, dans lequel il fonctionne selon une 
&equence choisie. et un mode en temps i/irtueL dans lequel il fonctionne de fagon 
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